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Abstract. We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various 
secondary particles (in addition to ;r'''s) with y's as the final decay state, that are produced in cosmic-ray (p's and a's) 
interactions with the interstellar medium. We derive an easy-to-use y-ray production matrix for cosmic rays with energies up 
to about 10 PeV. This y-ray production matrix is applied to the GeV excess in diffuse Galactic y-rays observed by EGRET, 
and we conclude the non-;r" decay components are insufficient to explain the GeV excess, although they have contributed 
a different spectrum from the ;r^-decay component. We also test the hypothesis that the TeV-band y-ray emission of the 
shell-type SNR RX J1713. 7-3946 observed with HESS is caused by hadronic cosmic rays which are accelerated by a cosmic- 
ray modified shock. By the statistics, we find a continuously softening spectrum is strongly preferred, in contrast to 
expectations. A hardening spectrum has about 1% probability to explain the HESS data, but then only if a hard cutoff at 
50-100 TeV is imposed on the particle spectrum. 
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y-ray production matrix in hadronic Interactions 

This work presents a careful study of the y-ray production in cosmic-ray (both p and a) interactions, by accounting 
for all decay processes including the direct production. For that purpose we employ the event generator DPMJET- 
ni lit] to simulate secondary productions in both p-generated and He-generated interactions. We include all relevant 
secondary particles with y-rays as the final decay products. For the composition of the ISM, we assume 90% protons, 
10% helium nuclei, 0.02% carbon, and 0.04% oxygen. Around the energy of production threshold, where DPMJET 
appears unreliable, we apply a parametric model |i2J, that includes the resonance production for the n production. We 
thus derive a y-ray production matrix for cosmic rays with energies up to about 10 PeV that can be easily used to 
interpret the spectra of cosmic y-ray sources. 

We consider all the decay channels and their decay fractions published by the Particle Data Group to account for 
all secondaries (resonances included) calculated by DPMJET and the parametric method. The calculation shows the 
non-Tr" resources in hadronic interactions have contributed about 20% of the total y-ray photons, mostly from directly 
produced y-rays and decays of Tj, Kf^ and |3]. 

In the cosmic -ray interactions, we calculate the y-ray spectrum contributed by decays of unstable secondary particles 

Qy{Ey) = Y,nisM I dEcR Ncr{Ecr) c^cr C7(£'c«) {Ej^Ecr) (1) 
k JEcR "^y 

where '^^^ is the y-ray decay spectrum from secondary species k. Eq. ([T]) can be re-written into 

Qy{Ei) =Y,nisMAEj NcR{Ej)c^ja{Ej) ^ ^{Ei,Ej) = Y."'SMAEjNcR{Ej)cpj ajM^j (2) 

k y j 

thus reducing this problem to a matrix operation with the y-ray production matrix Mfj for which, each element M,y 
shows the value of the resultant particle energy spectrum ^\Ey=Ei,EcR=Ej, with j being the index for the generating 
cosmic-ray particle (p or a) and / being the index indicating the y-ray energy. The energy binnings £, and Ej are 
defined with good resolutions IH]. 

We use the y-ray production matrix to analyze the observed spectra of diffuse Galactic emission and of the shell-type 
SNR RX J 17 13.7-3946. We find that the GeV excess is probably not the result of an inappropriate model of hadronic y- 
ray production. We also test the hypothesis that the TeV-band y-ray emission of SNR RX J1713. 7-3946 observed with 



HESS is caused by hadronic cosmic rays that have a spectrum according to current theories of cosmic-ray modified 
shock acceleration. 



Application: the GeV-band y-ray spectrum and the TeV-band spectrum of RX J1713.7-3946 

With the 7-ray production matrix, we calculate the diffuse 7-ray spectrum generated by the observed cosmic-ray 
spectrum |4]. Fig.[T](Left) shows the observed GeV-band 7-ray emission from the inner Galaxy ^ in comparison with 
the contributions from 71" decay as well as bremsstrahlung emission describe by a power-law spectrum <I>b(£'): 



with power-law spectral index Fg = 2.1, the electron energy density (Og = 0.1, 0.4, 0.8 eV/cm^, and the gas column 
density N[sm = 3 x 10^^, 8 x 10^^ 3 x lO^'cm^^^ respectively. The generating cosmic -rays are assumed with an 
energy density Pe = 0.75 eV/cm^. Models based on the locally observed cosmic-ray spectra generally predict a softer 
spectrum for the leptonic components, even after accounting for inverse Compton emission ||5|], so we may in fact 
overestimate the GeV-band intensity of the leptonic contribution. Nevertheless, it is clearly seen in this figure, that 
in the total intensity an over-shooting around Ey ~ 300-600 MeV appears in the modelled 7-ray energy distribution, 
whereas a deficit is present above 1 GeV The observed spectrum of diffuse emission is always harder than the model 
spectrum, and we therefore conclude that an inaccurate description of hadronic 7-rays is ruled out as the origin of the 
GeV excess. 

For the TeV-band 7-ray spectrum of the shell-type SNR J1713. 7-3946 observed by the HESS collaboration l'6'l, we 
use the 7-ray production matrix to test cosmic -ray acceleration models [7, 8], which predict a continuous hardening of 
the cosmic-ray spectrum up to a high-energy cutoff. We therefore parametrize the spectrum of accelerated hadrons as 



where is the step function and £0 = 15 TeV is a normalization chosen to render variations in the power-law index 
s statistically independent from the choice of spectral curvature, a. The cutoff energy, /Smax, is a free parameter. The 
normalization A^o is obtained by normalizing both the data and the model to the value at 0.97 TeV. By the x'^ statistics, 
we obtain the best-fitting values and the confidence ranges of the three parameters, Emax, and (7, given values as 
5 = 2.13, (7 = -0.25, and 

Emax ^ 200 TeV, i.e. no cutoff. The best fit, shown in Fig.[T](Right), involves a continuous 
softening and is thus not commensurate with expectations based on acceleration at a cosmic -ray modified shock [7, .8]. 
Noted that very valuable would be data in the energy range between 1 GeV and 200 GeV that may be provided by 
GLAST in the near future. Fig. |2] shows the confidence ranges of the parameters in Eq. (|4]l, with confidence levels of 
1, 2 and 3 sigma. The analysis strongly suggests a negative spectral curvature, (7 < 0, with confidence more than 95%, 
in contrast to the expectation of standard cosmic -ray modified shock models. 



We have considered a full picture of the hadronic 7-rays in cosmic -ray interactions and introduced an easy-to-use 7- 
ray production matrix which can be used for arbitrary cosmic-ray spectrum. The matrices are available for download at 
website http://cherenkov.physics.iastate.edu/gamma-prod We apply the production matrix to calculate the 7-ray GeV 
excess and also the TeV-band spectrum of SNR RX J1713. 7-3946. We conclude that 1) the modifications in the GeV- 
band 7-ray emission of hadronic origin are insufficient to explain the GeV excess in diffuse galactic 7-rays; 2) a soft 
cut-off at about 100 TeV is statistically required in the particle spectrum if the TeV-band spectrum of RX J 17 13. 7-3946 
as observed with HESS is caused by cosmic-ray nucleons; 3) no evidence for efficient nucleon acceleration to energies 
near the knee in the cosmic-ray spectrum, nor evidence of the spectral curvature and hardness predicted by standard 
models of cosmic -ray modified shock acceleration. We emphasize the need for GLAST data to better constrain the 
7-ray spectrum below 100 GeV. 
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FIGURE 1. Left: Diffuse y-ray spectrum at GeV range, shown in comparison witii data for the inner part of Galaxy at 315^ < / < 
345° and \b\ < 5° H. See text for discussions. Right: The TeV-band y-ray spectrum observed from RX J1713.7-3946 with HESS 
0], shown in comparison with the best-fit model of hadronic y-ray production of Eq. l[4j. 
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FIGURE 2. Left: The confidence regions for the spectral curvature a and the spectral index s. The shaded areas correspond to 
probabilities 68% (1 sigma), 95% (2 sigma) and 99.7% (3 sigma). Right: The confidence regions for the spectral curvature CT and 
the spectral index s, versus the cutoff energy iimax- The contour levels are arranged as the same. 
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